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Abstract 

In this letter, we investigate the magnetic properties, electronic structures, 
Slater-Pauling behaviours of some quanternary Heusler alloys with 4d and 
3d transition metal elements. The energy levels of the minority-spin elec¬ 
tronic band structures for all our calculated alloys are discussed. We have 
summaried 21 quanternary Heusler alloys. Our calculations indicate that the 
half metals ZrGoGrSi, ZrGoGrIn, ZrGoFeSi, ZrFeVAl, ZrFeVGa and NbFe- 
GrAl have large spin-flip gaps. 
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1. Introduction 

Efficient spin injection from a ferromagnet to a semiconductor is very 
meaningful for the development of the performance of spintronic devices 
[ij. Besides, since the prediction of half-metallic ferromagnetism (HMF) 
of Heusler alloy NiMnSb by hrst-principle calculations in 1983 j^, HMF has 
attracted great interest. And the HMF Heusler alloys are good candidates 
for the applications of spintronic devices, as these alloys often have high spin 
polarization, high Gurie temperature and compatible lattice structure. There 
are many investigations on the half metallic Heusler alloys. Lots of interesting 
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properties of them have been found . For example, the so called LiMgPdSn 
or Y-type structure [3, [sj Heusler compounds have attracted much attention 
for their potentional value. This structure Heuslers can be half metals (HM) 
and spin-gapless semiconductors (SGS). The quanternary Heusler ZrFeTiZ 
(Z=A1, Si and Ge) and ZrNiTiAl have been reported to be HMF recently 
with large band gaps and spin-flip gaps. 

Inspired by the above, we study a set of 21 quanternary Heusler alloys 
with 4d and 3d transition metal elements by using density function theory 
(DFT) calculations. The plan of our letter is as follows: The next Section 
2 gives the computational details; Section 3 presents the discussions of the 
calculated results; Finally the conclusions are given in Section 4. 

2. Computational details 

The lattice optimization, density of states, magnetic properties and band 
structure of our calculated quanternary Heusler alloys are calculated by DFT 
calculations. All our DFT calculations are performed by using the full- 
potential local-orbital minimum-basis band structure scheme (FPLO) [3,1^ 
with generalized gradient approximation (GGA)j^, [l^, [ll|. For the irre¬ 
ducible Brillouin zone, we use the k meshes of 20x20x20 for all the calcula¬ 
tions. The convergence criteria of self-consistent iterations is set to 10“® to 
the density and 10“^ Hartree to the total energy per formula unit. 

3. Results and discussions 

In our calculations, Xi is a 4d transition metal element Zr or Nb. X 2 and 
Y are the 3d transition metal elements Go, Gr, Fe, V and Ti, which the atomic 
number of X 2 is bigger than that of Y. Z is a main group element. LiMgPdSn 
or Y-type structure (space group No.216, F43m) have four Wyckoff positions 
4a(0, 0, 0), 4b(i,i,i), and 4d(|,|,|). Yi, Y 2 , Y and Z atoms can 

sit at three possible nonequivalent superstructures: Y-type (/) (4a, 4b, 4c, 
4d), Y-type (//) (4a, 4c, 4b, 4d) and Y-type {III) (4a, 4d, 4b, 4c). In order 
to get the equilibrium structures of the alloys, the geometry optimization 
is performed in their three different conhgurations by calculating the total 
energy as a function of lattice constant. According to the calculated results 
of total energy at equilibrium lattice constant, we can get that Y-type (/) 
in spin-polarization is the most stable structure for all our calculated alloys. 
So only the Y-type (/) structure in spin-polarization will be discussed. 
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The optimizated lattice constants, band gaps, spin-flip gaps and physical 
natnre are shown in Table 1. As can be seen, there are six SGSs. The detail 
investigation of these SGSs is in onr rprevions paper The discnssions 
of ZrGoFeSi and ZrGoFeGe are in onr another paper [l3|. As a contrast we 
have also stndied ZrFeTiAl, and onr resnlts £t well with the discnssions in 
Ref. j^. In onr calculations, there are 14 HMs in total. The 14 HMs have 
large band gaps. 

In Figure 1, the total spin magnetic moment {Mt/ ijlb) versus the number 
of total valence electrons in the unit cell is plotted. For the case of positive 
Mt, the spin-down and spin-up are the majority and minority states, while in 
the the case of negative Mt the opposite occurs. Our results can be classhed 
along two lines representing two variants of the Slater-Pauling (SP) rules. It 
can be seen the alloys consisted of Nb only comply with the Mt = Z* — 24 
SP rule. When Xi is Zr and Z is Be, Al, Ga or In, the quanternary Heusler 
alloys comply with the Mt = Zt — 18 SP rule. As for the case that Xi is Zr 
and Z is Si or Ge, the quanternary Heusler alloys are divided into two parts. 
One part obeys the Mt = Z* — 18 SP rule and the other complies with the 
Mt = Zt — 24 SP rule. Attention should be paid to that Zr and Ti are in 
the same Group IV B, while Nb and V are in the same Group V B. So the 
results of Zr is similar to that of Ti 
that of V 


14 , while the results of Nb is similar to 


To understand the SP behaviors of these quanternary alloys, we present 
a schematic representation of the energy levels of the minority-spin band 
structure for them in Figure 2. The occupied and unoccupied states lo¬ 
cate below and above the Fermi level, respectively. The symbol beside each 
orbital denotes the corresponding degeneracy while we have introduced the 
nomenclature of Ref. | jl5| . The origin of the SP rules for the Heusler alloys 
with 4d and 3d transition metal elements is due to the hybridization effects 
which have been discussed detailedly in our recent paper [l^. The Zr-based 
alloys which obey the Mt = Zt — 18 SP rule are similar to the Sc 2 - and Ti 2 - 
based inverse Heusler alloys jl^. As can be seen in Figure 2, the triple¬ 
degeneration tiu states and the double-degeneration e„ states are in such high 
energy level that both states are empty. And the minority-spin gaps are cre¬ 
ated between the non-bonding tiu and bonding t 2 g states. In the schematic, 
one single state Ixs, two triple-degeneration states 3xp as well as 3xf2g and 
one double-degeneration state 2 x are below the Fermi level. So the number 
of states below the Fermi level: Nj,=9. We can directly deduce the number 
of occupied majority-spin states: N j'= Ztot — N 1. Hence the total magnetic 
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moment: Mtot = (A^ t i)f^B = {Ztot - 2 x N Dfis = {Ztot - 18)/is. 
Therefore some Zr-based alloys obey the Mt = Zt — 18 SP rule. 

The case of some Zr- and Nb- based alloys which obey the Mt = Zt — 24 


SP rule is similar to the discussions of the full-Heusler alloys in Ref [15 


As shown in Figure 2, the tiu states are below the Fermi level, and the Cu 
states are still above the Fermi level. The tiu and states are occupied and 
empty respectively. As can be seen, the minority-spin gap is created between 
the nonbonding tiu and states. Below the Fermi level are the Ixs, 3xp, 
2xeg, 3xt2g and 3xti„ states. So NJ, is 12. Hence we can directly get that: 
Mtot = (A^ t Df'B = {Ztot - 2 X N Dub = {Ztot - 24 )hb. Therefore 
some Zr- and Nb- based alloys obey the Mt = Zt — 24 SP rule. 

In Figure 3, we show the total density of states (TDOS) for ZrCoCrSi, 


ZrCoCrIn, ZrCoFeSi [l3j, ZrFeVAl, ZrFeVGa and NbFeCrAl. For the other 
HMs, we can get similar conclusions. As can be seen, the TDOS in the 
majority-spin channel get through the Fermi level, which exhibits a metallic 
behaviour. And there is a band gap around the Fermi level in the minority- 
spin channel, which exhibits a semiconductive behaviour. Around the Fermi 
level, the energy gap in the minority-spin channel and the crossing over in 
the majority-spin channel lead to the 100% spin polarization. So these six 
alloys are HMs. The spin-flip gap is dehned as the minimum gap between 
the Fermi level and the edge of the band gap in the minority-spin channel. 
Therefore we can get that the spin-flip gaps are 0.10, 0.22, 0.22 13|, 0.33, 0.38 
and 0.21 eV for ZrCoCrSi, ZrCoCrIn, ZrCoFeSi [], ZrFeVAl, ZrFeVGa and 
NbFeGrAl, respectively. Hense they may keep their half-metallicity at room 
temperature. And they may be good candidates for the future spintronics 
use. 

We also discuss the partial density of states (PDOS) of ZrFeVGa as a 
representative. Figure 4 is the PDOS and TDOS. As can be seen, the 3d 
states of Fe and V make the main contribution to the TDOS in the range of 
-4 to 2 eV. The 4d states make the most contribution among all the states 
of Zr in this range. And the 4p states of Ga make the most contribution 
among all the states of Ga near the Fermi level. It is clear that there is a 
hybridization between the 3d states of Fe and V around the Fermi level. 


4. Conclusions 

In a word, we study a series of quanternary Heusler alloys with 4d and 3d 
transition metal elements in this letter. Specihcly, we study the electronic 
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structures, magnetic properties and Slater-Pauling behaviours of 21 alloys. 
The Slater-Pauling behaviours of these alloys have been summaried. The 
energy levels of them are studied in detail. We focus on six HMs which 
have large enough spin-flip gaps. The density of states of these HMs are 
investigated. These HMs may be stable at room temperature because of 
their large spin-flip gaps. And they may be good candidates for the future 
spintronics use. 
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Table 1: 

X 1 X 2 YZ 

^opt (-^) 

E,(eV) 


Nature 

ZrGoGrSi 

6.033 

0.27 

0.10 

HM 

ZrGoGrAl 

6.251 

0.93 

0.02 

HM 

ZrGoGrBe[12] 

6.013 

0.71 

0.26 

SGS 

ZrGoGrGa 

6.236 

0.89 

0.06 

HM 

ZrGoGrIn 

6.455 

0.83 

0.22 

HM 

ZrGoFeSi[l^ 

5.971 

0.64 

0.22 

HM 

ZrGoFeGe[13] 

6.056 

- 

- 

Nearly HM 

ZrGoFeP[12] 

5.944 

0.41 

0.06 

SGS 

ZrFeGrSi 

6.039 

0.52 

0.26 

HM 

ZrFeGrIn[12] 

6.419 

0.80 

0.06 

SGS 

ZrFeGrGa[12] 

6.184 

; o . 7 i t 0.02 

SGS 

ZrFeVAl 

6.248 

0.93 

0.33 

HM 

ZrFeVGe[12] 

6.210 

;0.81 t0.41 

SGS 

ZrFeVGa 

6.237 

0.98 

0.38 

HM 

ZrGoVIn[12] 

6.468 

0.98 

0.15 

SGS 

ZrFeTiAl[^ 

6.31 

0.56 

- 

HM 

ZrFeTiAl 

6.325 

0.80 

0.22 

HM 

NbFeVSi 

6.000 

- 

- 

Nearly HM 

NbFeVGe 

6.086 

0.18 

0.03 

HM 

NbFeGrSi 

5.912 

0.08 

0 

HM 

NbFeGrAl 

6.035 

0.43 

0.21 

HM 

NbGoGrAl 

6.026 

0.17 

0.06 

HM 
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Table Caption 

Table 1; The results of the lattice optimization {^opt), band gap (E^), spin- 
flip gap {Esfg) and Physical nature of our calculated quanternary Heusler 
alloys. 

Figure Captions 

Figure 1: Calculated total magnetic moment (Mj) as a function of the total 
number of valence electrons {Zf), in the unit cell for the calculated Heusler 
alloys. 

Figure 2; Schematic representation of the energy levels of the minority-spin 
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